Brain activity mapping from MEG data via a hierarchical Bayesian algorithm with automatic depth weighting: sensitivity and specificity analysis
نویسندگان
چکیده
A recently proposed IAS MEG inverse solver algorithm, based on the coupling of a hierarchical Bayesian model with computationally efficient Krylov subspace linear solver, has been shown to perform well for both superficial and deep brain sources. However, a systematic study of its sensitivity and specificity as a function of the activity location is still missing. We propose novel statistical protocols to quantify the performance of MEG inverse solvers, focusing in particular on their sensitivity and specificity in identifying active brain regions. We use these protocols for a systematic study of the sensitivity and specificity of the IAS MEG inverse solver, comparing the performance with three standard inversion methods, wMNE, dSPM, and sLORETA. To avoid the bias of anecdotal tests towards a particular algorithm, the proposed protocols are Monte Carlo sampling based, generating an ensemble of activity patches in each brain region identified in a given atlas. The sensitivity is measured by how much, on average, the reconstructed activity is concentrated in the brain region of the simulated active patch. The specificity analysis is based on Bayes factors, interpreting the estimated current activity as data for testing the hypothesis that the active brain region is correctly identified, vs. the hypothesis of any erroneous attribution. The methodology allows the presence of a single or several simultaneous activity regions, without assuming the knowledge of the number of active regions. The testing protocols suggest that the IAS solver performs well in terms of sensitivity and specificity both with cortical and subcortical activity estimation.
منابع مشابه
Hierarchical Bayesian estimates of distributed MEG sources: theoretical aspects and comparison of variational and MCMC methods.
Magnetoencephalography (MEG) provides millisecond-scale temporal resolution for noninvasive mapping of human brain functions, but the problem of reconstructing the underlying source currents from the extracranial data has no unique solution. Several distributed source estimation methods based on different prior assumptions have been suggested for the resolution of this inverse problem. Recently...
متن کاملAutomatic fMRI-guided MEG multidipole localization for visual responses.
Previously, we introduced the use of individual cortical location and orientation constraints in the spatiotemporal Bayesian dipole analysis setting proposed by Jun et al. ([2005]; Neuroimage 28:84-98). However, the model's performance was limited by slow convergence and multimodality of the numerically estimated posterior distribution. In this paper, we present an intuitive way to exploit func...
متن کاملEvaluation of hierarchical Bayesian method through retinotopic brain activities reconstruction from fMRI and MEG signals
A hierarchical Bayesian method estimated current sources from MEG data, incorporating an fMRI constraint as a hierarchical prior whose strength is controlled by hyperparameters. A previous study [Sato, M., Yoshioka, T., Kajihara, S., Toyama, K., Goda, N., Doya, K., Kawato, M., 2004. Hierarchical Bayesian estimation for MEG inverse problem. Neuroimage 23, 806-826] demonstrated that fMRI informat...
متن کاملHierarchical Bayesian estimation for MEG inverse problem.
Source current estimation from MEG measurement is an ill-posed problem that requires prior assumptions about brain activity and an efficient estimation algorithm. In this article, we propose a new hierarchical Bayesian method introducing a hierarchical prior that can effectively incorporate both structural and functional MRI data. In our method, the variance of the source current at each source...
متن کاملDynamical MEG source modeling with multi-target Bayesian filtering.
We present a Bayesian filtering approach for automatic estimation of dynamical source models from magnetoencephalographic data. We apply multi-target Bayesian filtering and the theory of Random Finite Sets in an algorithm that recovers the life times, locations and strengths of a set of dipolar sources. The reconstructed dipoles are clustered in time and space to associate them with sources. We...
متن کامل